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lO . ABSTRACT 

m ■ 

0^ Context. To constrain different cosmological models and especially the equation of state parameter w of dark energy, the phases and amplitudes 
of baryonic acoustic oscillations in the galaxy power spectrum can be used as cosmological probes. 

Aims. It will be shown that the phases as well as the amplitudes of baryonic acoustic oscillations can be extracted out of a galaxy power spec- 
' trum, including various observational effects like growth suppression, redshift space distortions, and galaxy biasing. The phases of baryonic 
, acoustic oscillations are used as a standard ruler for a cosmological test to constrain w over a large redshift range with a minimum of assump- 
^ l' tions. 

Methods. A non-oscillating phenomenological fitting function is used to extract the oscillatory part of the galaxy power spectrum and to dis- 
entangle phase information from amplitude information. The method (FITEX) is tested with simulated data of the Hubble Volume Simulation 
^ to include redshift space effects, non-linear structure growth and biasing. A cosmological test is introduced, which compares the extracted 

oscillations against a theoretical model template to derive constraints on the w parameter. 
|L* Results. The phenomenological fitting function is found to model the various distortions of the galaxy power spectrum to the sub-percent level. 
. The various distortions only boost the amplitude of the oscillations. The theoretical template is accurate enough to test for small deviations in 
the phases of the oscillations, resulting from different w values. A cosmological test, using the baryonic acoustic oscillations as a standard ruler, 
5_j is able to constrain w in a robust way. These constraints are based solely on the phase of these oscillations, assuming no errors of the value 
at the end of the Compton drag epoch, as the standard ruler. The w constraints are expected to further improve when the full two dimensional 
redshift space effects and the amplitudes of the baryonic acoustic oscillations are taken into account. 
Conclusions. 
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1. Introduction 

Acoustic oscillations as observed in the temperature 
anisotropics of the cosmic microwave background radia- 
tion (CMB) are traditionally used to constrain the values of 
certain cosmological parameters. The intrinsic amplitudes and 
locations of the oscillations are determined by the densities and 
pressures of the various energy components in the very early 
and hot Universe. The oscillations are furthermore modified by 
their subsequent geometric projection onto the present hyper- 
sphere, where they are observed. It turns out that both intrinsic 
CMB oscillations and their geometric projection are mainly 
shaped by the physical properties at high redshifts. Excluding 
the very large scales where cosmic variance makes analysis 
quite difficult, the observed oscillations in the CMB are highly 
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degenerate against certain changes of the energy density poE 
and the equation of state parameter, w - Pde/Pdec^, of the 
dark energy, dominating at low redshifts (see also Caldwell, 
Dave & Steinhardt 1998). 

However, related so-called baryonic acoustic oscillations 
(BAOs) are expected to be observable at lower redshifts in the 
matter power spectrum. Therefore, the phases of the BAOs can 
be used as cosmic rulers in a similar manner as the CMB os- 
cillations, but now at much smaller redshifts (Eisenstein, Hu 
& Tegmark 1998). Here, the observed phases are affected by 
the late-time geometry and thus by the value of w, the prime 
cosmological parameter of the present investigation. 

The BAOs are classical Doppler peaks in the density dis- 
tribution of matter (e.g. Hu & Sugiyama 1996). They are trig- 
gered in the hot Universe by oscillatory velocity patterns of 
the baryons. These velocity oscillations are generated, in the 
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same environment as the CMB density oscillations, by sound 
waves on scales where radiation pressure could stabilize the 
fluid against gravitational collapse. At later times, this oscil- 
latory velocity pattern of the baryon field produces kinemati- 
cally a new field of matter density fluctuations in form of quasi- 
regular matter oscillations. The amplitudes of these BAOs are 
determined by the ratio of baryonic matter to the overall matter 
density. The BAOs thus constitute a quasi-regular pattern of os- 
cillating substructures superposed with small amplitudes onto 
a general density field which fluctuates irregularly with several 
100 times higher amplitudes. 

The imprints of the BAOs on galaxy and clusters distribu- 
tions, gravitational shear maps etc. can provide at least in prin- 
ciple a clean cosmic ruler for precise tests of the w parame- 
ter The characteristic scale, s, of these imprints is the comov- 
ing distance, sound waves can travel during the epoch where 
baryons and photons are strongly coupled through Compton 
and Thomson scattering (Compton drag). Cosmological tests 
based on the resulting sound horizon, s, as a metric ruler are 
expected to have only very small systematic errors because the 
phases of the BAOs located at scales small compared to s are 
solely determined by well-understood physical processes. A 
practical problem is, however, to separate the oscillatory BAO 
modes from the irregular fluctuation field with sufficient accu- 
racy. 

One might object that the small amplitudes of the BAOs can 
easily be washed out especially by structure growth which can 
mix perturbation modes with adjacent wavenumbers k. This is 
certainly true for oscillations in the non-linear regime (Meiksin, 
White & Peacock 1999). On larger scales, however, recent ob- 
servations of 2dF galaxies at redshifts z < 0.3 clearly show 
several BAOs in the galaxy power spectrum (Cole et al. 2005). 
In addition, the sequence of BAOs in A:-space is projected into 
a single wiggle in the space domain. The corresponding ex- 
cess galaxy correlation was in fact observed in the two-point 
spatial correlation function of the luminous red SDSS galax- 
ies at z < 0.5 (Eisenstein et al. 2005). In addition to these 
observational indications, results from recent numerical sim- 
ulations suggest that within certain redshift and scale ranges, 
stable BAOs could exist as useful probes for cosmological in- 
vestigations (e.g. Seo & Eisenstein 2005, Springel et al. 2005, 
Jeong & Komatsu 2006). However, several critical issues like 
the exact behavior of baryons during structure growth on BAO- 
scales can only be discussed with much larger simulations and 
a better understanding of galaxy formation. 

Applications of BAOs for cosmological tests of w are con- 
fronted with the following situation. Present observational con- 
straints of w using different combinations of CMB, galaxy clus- 
ters, galaxies, and gravitational lensing data are all found to be 
consistent with w - -I, i.e., the value for the cosmological 
constant. The Icr error of the w values, derived from partially 
dependent observations, is around 10-20% (for a recent review 
see, e.g., Schuecker 2005). To improve current estimates, we 
thus have to extract the BAOs and to measure their phases quite 
accurately in the presence of non-linear and scale-dependent 
effects. Therefore, even tiny systematic errors in the analysis 
on percent levels can have severe consequences for the obser- 
vational accuracy of w. 



The basic aim of the present paper is to describe a new, sim- 
ple, and robust method to extract BAOs from (non)linear, scale 
and redshift-dependent biased, redshift-space galaxy power 
spectra under realistic lightcone survey conditions. For brevity, 
we call the method 'fit and extract' (FITEX). The method has 
the potential to reduce many of the abovementioned sources of 
systematic errors. Our starting point is to use as few assump- 
tions as possible for the extraction of the BAOs from a complex 
power spectrum. For reasons which will become clear later, we 
do the separation by fitting a flexible non-oscillating function 
to an observable which is directly related to the transfer func- 
tion and not to the power spectrum itself. The crucial point is to 
show that BAOs extracted in this simple manner are still gov- 
erned by simple physical processes, and that FITEX stays ro- 
bust even under complex survey conditions. We believe that 
such model-independent approaches are quite important, espe- 
cially in light of the fact that at least in the near future we can- 
not expect to model the abovementioned non-linear and scale- 
dependent effects with accuracies on the sub-percent level. 

Our computations are in several cases optimized to galaxy 
surveys covering typically several 100 square degrees on the 
sky in the redshift range between z = 2-4. This choice 
is motivated by the Hobby Eberly Dark Energy Experiment 
(HETDEX), which is planned to measure the w parameter with 
several million Ly-o' emitting galaxies at these redshifts (Hill 
& McQueen. 2004, Gebhardt et al. in prep.). We see the results 
of the present paper as a useful contribution to a realistic error 
forecast for this important project. 

We organized the paper as follows. In Sect.|2lwe motivate 
the basic idea of the FITEX algorithm. In the following sec- 
tions we test its performance to separate BAOs from a complex 
power spectrum under linear conditions (Sect.|3}, under quasi- 
non-linear conditions (Sect.0}, in redshift space (Sect.|5}, and 
from biased samples (Sect.|6}. Finally, we illustrate the perfor- 
mance of FITEX including all eff'ects and use a cosmological 
test of w as a benchmark to investigate the quality of the BAO 
extraction together with the theoretical template (Sect.0. Our 
results are mainly based on a deep data wedge extracted by the 
Virgo Consortium from the Hubble Volume Simulation. The 
paper thus presents for the first time constraints on the appli- 
cation of BAOs for cosmological tests under realistic lightcone 
survey conditions. 

If not mentioned explicitly, we assume a spatially flat 
Friedmann-Lemaitre Robertson- Walker world model with the 
Hubble constant in units of h = //o/(100kms"' Mpc"'), the 
present values of the total matter density Q^h^ = 0.147 and 
baryon density Q^h^ = 0.0196, the density of relativistic mat- 
ter (e.g. neutrinos) O,, = 0, and the mean CMB temperature 
T^CMB - 2.728 K (concordance cosmology). 

2. Baryonic Acoustic Oscillations (BAOs) 

In this section, we discuss the 'theoretical wiggle function', 
that is, a reference function (see Eq.[n below) we use in our 
cosmological tests of w to match the BAOs extracted in a 
certain manner from the complex galaxy power spectrum. 
The power spectrum may be written in the form P{k,z) - 
Aik, z, b) k" T^{k). The amplitude A{k, z, b) includes the primor- 
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dial amplitude, redshift and scale-dependent effects of redshift 
space distortions, linear and non-linear structure growth, and 
galaxy biasing which will be discussed in the course of the 
paper The exponent n is the slope of the primordial power 
spectrum, and T{k) the transfer function. Our further treatment 
of the formation of BAOs is based on Hu & Sugiyama (1996) 
where a detailed description of the relevant physical processes 
can be found, and on the fitting equations derived by Eisenstein 
&Hu(1998). 

Between the epoch of matter-radiation equality at redshift 
Zeq ~ 3526 (concordance cosmology) and the end of the 
Compton drag epoch at Zd ~ 1026, the baryons follow a reg- 
ular velocity pattern which can generate new density fluctua- 
tions kinematically. On small scales and at Zd, this effect (veloc- 
ity overshoot) overrides the intrinsic density fluctuations of the 
baryons. As adiabatic modes dominate the isocurvature modes, 
the baryon density oscillates as sin(A:s). The phase of the BAOs 
is frozen out at Zd at the value ks, with 5 « 152 Mpc being the 
comoving sound horizon at Zd and k the comoving wavenum- 
ber of the oscillation. The BAOs are thus n/l out of phase with 
the corresponding CMB fluctuations. The sound horizon s at Zd 
is the standard ruler we are looking for Its value calibrates the 
phases of the theoretical wiggle function and can be measured 
quite easily with CMB experiments on the sub-percent level. 
This includes the contribution of unknown relativistic energy 
components (Eisenstein & White 2004). For consistency with 
the Hubble Volume Simulations (see Sect. ^ we use Ov = 0. 

The sinusoidal fluctuations in the baryon density are damp- 
ened by expansion drag, gravitational forcing, and Silk damp- 
ing. In addition, cosmic expansion forces the velocity contri- 
butions to fall off at large scales, and amplitudes decline when 
dark matter dominates the energy density. While these effects 
only reduce the amplitude of the fluctuations and do not af- 
fect the phase, the damping may be summarized by the pseudo 
transfer function T„(k) ~ ;o(/ts)e"<*^*^'"')"7[l + (fiblksf], 
with Jo the spherical Bessel function of order zero, the ex- 
ponent nis ~ 1.4 which is basically independent from cos- 
mology, ^siik = 1.6(Qb/iy (Om /2')°"[l + (10.4Q „/2^)-°''^], 
and Pb = 0.5 H- /b -H (3 - 2/b) ^J{V^2nJ^^^ with /b = 
Q.bl^m- However, velocity overshoot dominates only on scales 
small compared to s. On larger scales, the original sound hori- 
zon, s, appears to be reduced by {sj sf = 1-1- (J3/ks)^, with 
P = 8.41(Q,n/!2)0-435. This con-ection of s is about 0.2% at 
60 /z"' Mpc and 1 .3% at 120 /z ' Mpc. Though relatively small, 
the corrections are important as they change the length of the 
standard ruler s depending on k. Collecting all scale-dependent 
terms we get the oscillatory solution (theoretical wiggle func- 
tion) 



pressureless and will fall into the potential wells of CDM. This 
results in a transfer function valid for both species of matter. 



TAk) 



1 + (Pb/ks) 



3 Mk s) . 



(1) 



A way to extract BAOs from a complex power spectrum 
can be found when we specify the relation between the wiggle 
function, T„{k), and the total transfer function T(k), which we 
discuss now. Each particle species, in the present case CDM 
and baryonic matter with the corresponding densities and 
D.b, should have separate effective transfer functions, and 
Tb- Though, after the drag epoch at Zd, baryons appear basically 



Oc Qb 
T(k) = -^TAk) + -^Tbik), 



(2) 



with Tc{k) - fT()(k, 1, /?£.) + (! -f)To(k, acPc), written in terms 
of the generalized transfer function f()(k,ac,j3c) as defined in 
Eisenstein & Hu (1998). Here, / = 1/[1 -h (ks/5.4)'^] smooflis 
the combination of the almost baryon-free and baryon-loaded 
solutions near s, and 



Tbik) = fb(^) + TAk), 



(3) 



where f\,(k) 



^"i+iksfsi)^^ ■ These approximations are better 
than 2% for Qb/f^m < 0.5. Note that Ob/f^m corre- 
sponds to ffcySc — » 1. If Qc >> Qb, the effects of dark mat- 
ter dominate over velocity overshoot of the BAOs. Written in 
this way, we immediately see that the oscillatory part Q is 
added on top of a non-oscillatory part and can be computed by 
subtracting off a smooth continuum from the complete trans- 
fer function. This process establishes the basic methodology of 
FITEX. 

The different terms of the transfer function are shown in 
Fig-ID The upper panel shows the components which determine 
the global shape of the transfer function. Note that an additional 
oscillating component, the first term in Eq. Q, is necessary to 
describe the transition from large to small scales with sufficient 
accuracy. This is a direct consequence of the fact that the main 
effect of the baryons on the transfer function is a damping of 
the overall growth of dark matter between Zeq and Zd which sig- 
nificantly reduces the fluctuation power at scales smaller than s. 
The resulting break is the imprint of the baryons on the power 
spectrum which is easiest to observe. The term fb thus con- 
tributes mainly to the shape as it is rapidly declining on scales 
smaller than s to model the damping effect of the baryons. 

In the lower panel of Fig. [2 the two parts of the baryon 
transfer function are shown. One can clearly see the first to 
fifth oscillations. Note that Tb declines very fast and is almost 
zero after the second oscillation. It is compared to the second 
term of the baryonic transfer function, T^ik), and a reference 
function, r„yef, which is formally constructed to contain all os- 
cillations. Thus, the reference function includes both the small- 
scale oscillations from T^ik) and the large-scale oscillations 
from f\,{k). The latter is, however, also determined by the dark 
matter function T(i{k, 1,1) (see below Eq.|3|l. Therefore, exclud- 
ing or including knowledge of the shape of the dark matter 
function for the extraction of BAOs corresponds to working 
either with T^ik) or with r„,_ref ■ Both T^ik) and r„,_ref agree rea- 
sonably well beyond k > 0.05 /i Mpc ', especially when one 
keeps in mind that only the phases of the oscillations are rel- 
evant as the amplitudes are subject to distortions anyway (see 
below). 

The good agreement between the phases of T„{k) and r„_ref 
means, that quite simple physics as described by (mainly 
velocity overshoot and Silk damping) dominate the oscillations 
expected to be seen in the overall power spectrum. However, on 
larger scales, the physics get more complex due to the grow- 
ing influence of CDM and the backreaction of the baryons 
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3. Extracting BAOs under Linear Conditions 
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Fig. 1. Comparison of different terms of the transfer function com- 
puted with fib 'i^ = 0.0196, £l,^/i2 = 0.147 and TcMB = 2.728 K. 
Upper panel: Non-oscillatory parts. The function r^hape is the sum of 
the other three. Lower panel: Oscillating parts. The residuals of the 
subtraction between the transfer function and its corresponding non- 
oscillatory fit is shown as long-dashed line. Note the different scaling 
of the y-axes. See main text for more details. 



on structure growth. Thus r„, and r„, ief differ accordingly. 
Nevertheless, on small scales, r„, is a good theoretical wiggle 
function for BAOs phases which can be used as a standard ruler 
for cosmological tests, with a minimum of theoretical assump- 
tions. 



The main advantages of constraining cosmological parameters 
with BAOs are its small systematic errors, the ability to dis- 
criminate between geometrical effects (homogeneous universe) 
and structure growth effects (inhomogeneous universe) of cos- 
mological parameters, and the potential to constrain dark en- 
ergy without assuming a certain dark matter model. The ob- 
served phase of the BAOs is only affected by the geometry of 
the universe, while the observed amplitude of the oscillations is 
mainly affected by structure growth and other amplitude effects 
(see Sects. 0]to|6j. It is thus very important to extract the BAOs 
from the power spectrum in a way that does not mix these two 
benchmarks. As suggested by the previous discussion, BAOs, 
described by simple physical processes, should be extracted by 
subtracting the shape of the non-oscillatory part of the transfer 
function. This method is able to disentangle phase effects from 
amplitude effects in an effective manner 

In contrast to this method, previous studies (Blake & 
Glazebrook 2003, Angulo et al. 2005, Hutsi 2005, Seo & 
Eisenstein 2005, White 2005) used the CDM-based theoretical 
non-oscillatory power spectrum, Pretik, z), given in Eisenstein 
& Hu (1998) to analyze the BAOs. This is done by divid- 
ing the measured power spectrum, Pobs(^, z), by this reference 
power spectrum. Assuming that the transfer function can be 
split into a shape part, Tshape, and an oscillatory part, Twiggie, 
where Tshape > T^wiggie (Sect.|2j, the result of such a division is 



-Pobs(^,Z) Aobs(k,z)[Tshape(k) + ^wiggle 



7"wiggle (k) A,Uk,z) 

Ai-ef(^, Z) rshape(^) Ai-ef(fe, z) 



(4) 



Eq.Q gives a spectral ratio, which strongly deviates from the 
simple functional form of Eq. ([0 expected from basic physi- 
cal principles, and thus unnecessarily complicates any further 
analysis of the resulting "BAOs". A possible result of such a 
division is given in Springel et al. (2005, Fig. 6). It can easily 
be seen that phase information and amplitude information is 
mixed. 

In addition, the use of a theoretical non-oscillatoy reference 
power spectrum, Pi-efik, z), or the related theoretical "boosted" 
transfer function, A^gf(k, z)T^ef(k, z) = ■\/Prefik,z)/k", has a 
number of drawbacks: 

(1) The shape of the transfer function and thus the power 
spectrum vary with the cosmological parameters. To compute 
a reference power spectrum one has to know the exact values 
of several cosmological parameters. Furthermore one has to as- 
sume a CDM model as a prior (2) The reference power spec- 
trum has to be flexible to compensate for various distortions: 
structure growth, redshift space distortions, biasing etc. distort 
the transfer function and thus its shape. Most of the effects vary 
with redshift as well as with k. (3) The reference power spec- 
trum has to be exact on the sub-percent level. The BAOs make 
up only ~ 2% of the transfer function depending on the ra- 
tio Ob/^^m- An analytical function has to be more accurate to 
extract the oscillations from the transfer function. (4) It is non- 
trivial to provide an assumption-free template for the resulting 
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power ratio to perform a cosmological test. Instead, very accu- 
rate knowledge of the cosmology as well as amplitude effects 
is needed. 

Most of these issues are well-known and useful corrections 
are already available from the simulations. However, these cal- 
ibrations have certainly not reached the sub-percent accuracy 
level which is necessary for the cosmological tests, and one 
can doubt whether this is achievable at all. Furthermore, one 
of the main advantages of the BAO method, the ability not to 
assume a certain dark matter model, is negated. 

The Eqs. ( Ill3t instead suggest that the computation of the 
difference of transfer functions to extract BAOs is the most 
direct approach. Therefore, FITEX does not divide the ob- 
served power spectrum by a theoretical non-oscillating refer- 
ence power spectrum but subtracts a phenomenological non- 
oscillating fitting function from the observed boosted transfer 
function, to extract the BAOs. Formally FITDEX can be de- 
scribed as follows: 



J^^^-F(k,b,z) = yfAik:Zb)^TUk). (5) 
V k Qm 



We found that the formula 



F{k) = 



Jk/hf 



(6) 



is able to fit all non-oscillating distortions to a satisfactory de- 
gree and leaves enough free parameters to allow for a wide 
range of transfer functions. Eq. (|6jl has no oscillatory compo- 
nents and can thus only trace the shape of a transfer function. 
Most of the drawbacks of the reference power spectrum divi- 
sion can be avoided: 

(1) No cosmology or dark matter model has to be assumed 
to extract the BAOs. (2) The formula is flexible enough to com- 
pensate for various non-oscillating distortions without having 
to care about the source or physics of the distortions (see Sect. 
|3-|6|l. (3) The formula is able to fit power spectra reliably on 
the sub percent level (see Sect.|^. (4) It is possible to provide 
a theoretical template function that can be used without assum- 
ing knowledge about dark matter and various amplitude effects 
by leaving the amplitude as a free parameter, as phase informa- 
tion and amplitude information is disentangled. The amplitude 
could be used for a further cosmological test. 

To test FITEX, i.e., the combination of both the fitting func- 
tion F(k) and the wiggle function r„ (A;) under linear conditions, 
we first computed normalized transfer functions with CMBfast 
with the parameter values w = -1 (the transfer functions are 
affected by w only on scales above several Gpc, Ma et al. 
1999), Qm = 0.3, Qb = 0.04, h = 0.70, the primordial slope 
n = I, the average CMB temperature Tcmb = 2.728 K, the 
He mass fraction after primordial nucleosynthesis Ym - 0.24, 
and the number of neutrino families A^v = 3.04. In the second 
step, the phenomenological continuum function (Eq.|6|l was 
fitted to the multi-component transfer functions in the range 
O.OI < k < 0.3h Mpc~'. The results were subtracted from the 
original transfer functions to extract the BAOs. One example 
with f2m - 0.3 is shown in Fig.|2 In the first panel, the multi- 
component transfer function is plotted (solid line) as well as the 



phenomenological fitting function (dashed line). In the second 
and third panels, the differences between the transfer function 
and the continuum fit are plotted (solid lines), together with 
the theoretical wiggle functions (dashed lines) computed with 
Eq.lUli. 
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Fig. 2. Performance of the phenomenological fitting function under 
hnear conditions. The upper panel shows the CMB fast transfer func- 
tion (sohd line) and the best non-oscillating fit (dashed line). In the 
lower panels, the residuals of the subtraction are plotted (solid line) 
as well as the theoretical transfer function (dashed line). The original 
model gives the performance of the original theoretical function (Eq. 
Q while the corrected model uses the modified function (Eq.Q 

With typical baryon fractions, ft, ~ 0.15, baryonic oscil- 
lations are hardly visible as they make up only ~ 2% of the 
multi-component transfer function. This illustrates the main 
observational challenge in using wiggles as a metric ruler for 
cosmological investigations. On small scales, k > 0.15/;Mpc"', 
Fig-Elshows that the theoretical wiggle function describes the 
baryonic oscillations to a satisfactory degree. This is expected, 
as the transfer function is dominated by the effects of veloc- 
ity overshoot and Silk damping on these scales and the wiggle 
function was designed to model these effects. On large scales, 
k < 0.05 /iMpc ', the BAOs have to cope with competing CDM 
effects that dominate in these A;-ranges. As the influence of 
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CDM increases, the transfer function can not be described by 
baryonic physics alone and further corrections and assumptions 
have to be applied. Two related effects are: 

(1) Turnover-effect: The turnover in the matter power spec- 
trum occurs on the scale k^q of the particle horizon at Zeq, which 
coincides with the approximate location of the first wiggle. No 
observation has hitherto clearly revealed this turnover in spec- 
tral power The main problem in identifying the turnover is 
mostly due to the fact, that it is located at large scales where 
precise measurements of spectral power are difficult: the sur- 
veys cover too small volumes with slice-like shapes leading 
to strong smoothing (up to factors of two for the 2dF sur- 
vey, Percival et al. 2001) and significant leakage in the de- 
rived power spectra. Under such conditions it appears quite 
difficult if not impossible to discriminate the first wiggle from 
the turnover when even the position of the turnover could not 
yet been determined to some accuracy. To do this, an excel- 
lent signal to noise ratio would be needed as well as precise 
knowledge of the underlying linear theory matter power spec- 
trum. This includes information of distortions of the power 
spectrum to the sub-percent level. As a consequence, a model- 
independent approach like separating the BAOs from the multi- 
component transfer function with a phenomenological fitting 
function will only be able to detect the first wiggle when the 
baryon fraction is about as high as 40 percent or more, as 
we found with CMBfast simulations. However, we will show 
that our cosmological test is most sensitive to changes of w 
on scales around k = 0.1/iMpc Therefore, the inability of a 
model-independent approach to detect the first BAO does not 
matter much as this first wiggle barely contributes any infor- 
mation to the cosmological test, even more as the amplitude of 
the first wiggle is comparatively low as it is dampened by CDM 
effects (see Sect.0. 

(2) Phase-shift effects: As mentioned in Sect.|2 the in- 
creasing effect of dark matter on large scales introduces a k- 
dependent phase-shift of the BAOs. The theoretical wiggle 
function (Eq.^ already includes a phenomenological correc- 
tion for this effect (replacing s by s). A similar effect on the 
phases, but in the opposite direction, is introduced by the damp- 
ing of baryons on structure growth. In fact, the main effect of 
baryons on the multi-component power spectrum is, assum- 
ing a standard cosmology, a sharp break in power starting on 
scales smaller than the turnover. Formally, this break is caused 
by the rapid decay of the Ti, component of the transfer func- 
tion (see Fig. Q. As the break occurs right after the turnover, 
it is very hard to distinguish shape components from oscilla- 
tions. For the extraction of BAOs in that scale range we have 
found that this spectral break must be modeled in detail (with 
an accuracy on the sub-percent level) or artificial phase-shifts 
with sizes of the order of interesting w effects are introduced 
when naive fitting functions ignoring this baryonic backreac- 
tion are subtracted from the multi-component power spectrum. 
Thus the model independent fitting function presented in this 
work (Eq.|6jl, but also traditional wiggle-free dark matter trans- 
fer functions with a slightly different functional approach like 
those given in Bardeen et al. (1986) or Efstathiou, Bond & 
White (1992), introduce artificial phase-shifts of the generic 



form 



3/2 



0.98 -I- 



5.1 (Q,,/!^' 



2^0.47 



ks 



3/2 



(7) 



which gives a useful summary of the phase-shifts (about 0.7% 
at 60/!"'Mpc and 4% at 120/i"'Mpc) even for very different 
baryon fractions. This equation was derived phenomenologi- 
cally from fits to CMBfast simulations. The second panel in 
Fig-El shows ™ example for the concordance cosmology. It is 
seen that the phase-shifts are most prominent at intermediate 
scales, k < 0.15 /iMpc"', and counteracts the phase-shifts in- 
troduced by large-scale CDM potentials. The dashed line in 
the fourth panel shows the corrected theoretical wiggle func- 
tion, replacing s in Eq. Q by i from Eq. 0. It is clearly seen 
that the extracted wiggle function follows the predicted theo- 
retical, and phenomenologically corrected, wiggle function on 
the 0.5%-level in the range 0.05 < k < O.S/zMpc ' which is 
relevant for a cosmological test of the parameter w. In the fol- 
lowing, Eqs. Q and constitute our final theoretical wiggle 
function which we use in FITEX to match the BAOs extracted 
from the simulated data. 



4. Extracting BAOs under Quasi-Nonlinear 
Conditions 

The simulated data used to test FITEX for BAO extraction un- 
der more realistic conditions are provided by the ACDM ver- 
sion of the Hubble Volume Simulations conducted by the Virgo 
Consortium (Evrard et al. 2002). One billion dark matter par- 
ticles were simulated with h = 0.7, 0^ = 0.3, Qa = 0.7, 
f2b - 0.04, Qy = 0, and the normalization of the matter power 
spectrum, erg - 0.9, in a cube with a comoving length of 
L = 3000 h-^ Mpc and a mass of 2.25 ■ 10'^ h'^ Mq per par- 
ticle. The simulations were started with a 'glass-like' load (see 
e.g. Baugh et al. 1995). During the simulation itself, long range 
gravitational forces were computed on a 1024^ grid yielding a 
Nyquist critical frequency of kc - 1.07 /iMpc"'. The short 
range gravitational forces were computed via direct summa- 
tion and softened on a scale of 0.1 /i ' Mpc. We used the 
10x10 deg^ fraction of the XW extended deep wedge which 
uses periodic boundary conditions from a redshift of z = 4.4 to 
provide a survey lightcone up to the redshift Zmax - 6.8. The 
lightcone output comprises of data that includes cluster evo- 
lution and thus mimics real life observations. Larger distance 
from the observer correspond to higher redshifts where struc- 
tures are less pronounced due to linear and non-linear growth. 

Power spectra, Pobsik, z, b), were estimated at different red- 
shifts and for different galaxy biasing schemes b (see below) 
with the Fast Fourier Transform using the variance-optimized 
method of Feldman, Kaiser & Peacock (1994) for cubes along 
the line-of-sight (LOS) of the XW wedge. All the cubes were 
selected to fully fit into the lightcone in order to not intro- 
duce a non-trivial window function. The parameters of the 
cubes are given in Tab.[l] Note that the spectroscopic SDSS 
survey of normal galaxies covers an effective volume which 
is comparable to cubes with L ^ 500/;"' Mpc. The errors 
of the power spectral densities at wavenumber k were esti- 
mated with standard mode counting arguments by crplPik) - 
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Table 1. Parameters of the cubes for the power spectrum esti- 
mation. Col. 1 : Cube number. Col. 2: Comoving length of the 
cube (concordance cosmology). Col. 3: Redshift range covered 
by the cube. Col. 4: Fundamental mode of the discrete Fourier 
transformation which corresponds to the sample bin size of the 
power spectrum in k space. Col. 5: Nyquist critical wavenum- 
ber of the Fast Fourier Transformation. 



Sample 


length 


z-range 


Ak 


kc 




[h-' Mpc] 




[A Mpc"'] 


[AMpc"'] 


cubel 


200 


0.58 - 0.68 


0.031 


4.0 


cube2 


226 


0.68 - 0.79 


0.028 


3.6 


cube3 


256 


0.79 - 0.92 


0.025 


3.1 


cube4 


291 


0.92-1.1 


0.022 


2.8 


cube5 


329 


1.1 - 1.3 


0.019 


2.4 


cube6 


373 


1.3-1.6 


0.017 


2.2 


cube? 


423 


1.6-2.0 


0.015 


1.9 


cubeS 


479 


2.0-2.5 


0.013 


1.7 


cube9 


500 


2.5-3.2 


0.013 


1.6 


cube 10 


500 


3.2-4.1 


0.013 


1.6 


cube 11 


500 


4.1-5.4 


0.013 


1.6 



'\/2n{ 1 + 1 / [P(k)n] ]/kL to provide a first estimate. Here, n is the 
mean comoving particle number density in the box, and L the 
length of the box. We found that the initial glass load reduces 
the shotnoise level significantly below the standard Poisson 
case at redshifts z > 2 (see Smith et al. 2003). Instead of cor- 
recting our spectral estimator for this artifact (see, e.g.. Smith 
et al. 2003), we could in the end neglect its impact on the cos- 
mological test, because our phenomenological fitting function 
was flexible enough to compensate for the small residuals of 
this artifact seen in our standard scale range k < 0.3 li Mpc ' . 

In order to extract the BAOs from the power spectra 
Poh!,{k,z,b), we computed the corresponding boosted transfer 
functions first, ■\JA(k,z,b)T(k) = yjPohsik, z, b)/k, assuming 
a primordial spectrum with n - I (see below). Fig. shows 
the estimated power spectra and the corresponding boosted 
transfer functions for an unbiased dark matter particle distribu- 
tion in configuration space. Error bars include shot noise and 
sample variance. In most spectra the 2nd and 3rd BAOs are 
clearly visible. For k > 0.05 /i Mpc ', the phenomenological 
fitting function (Eq.|6|l can be computed very accurately as it is 
constrained by a large number of spectral data with relatively 
small errors, that can be traced reasonably well by a power law. 
Sample variance is negligible, as in this A:-range, each spec- 
tral point represents the average over a large number of modes. 
For k < 0.05 /i Mpc only a very small number of spectral 
densities, each estimated with a small number of independent 
modes, exists and the phenomenological function has to fit a 
wiggle that coincides with the turnover Therefore, the results 
in the range < A: < 0.05 h Mpc ' vary quite a bit. However, 
we already excluded this range from the cosmological test (see 
vertical dashed lines in Fig.|3land the discussion at the end of 
Sect.|3ll. 

The extracted BAOs for the four cubes 7- 10 of the XW deep 
wedge simulation corresponding to the redshifts z = 1.8-3.7 



L [h-' Mpc] 




0.01 0.05 0.1 

k [h-i Mpc]-' 



L [h-' Mpc] 
500 100 50 




0.01 0.05 0.1 

k [h-i Mpc]-i 



Fig. 3. Upper panel: Power spectra of the matter distribution from 
the cubes 7 (upper spectrum) to 10 (lower spectrum) of the XW deep 
wedge N-body simulation with superposed linear theory predictions. 
The linear matter power spectrum was taken from Eisenstein & Hu 
(1998). Amplitude corrections for linear structure growth were im- 
plemented in the standard manner. For better illustration, the 2nd, 
3rd and 4th spectrum is shifted downwards by the factors dex(O.l), 
dex(0.2), dex(0.3), respectively. The extraction of BAOs described in 
the main text is restricted to the A;-range bordered by the two vertical 
dashed lines. Lower panel: Same as upper panel for the correspond- 
ing boosted transfer functions with superposed fits of the phenomeno- 
logical fitting functions. The fit at z = 2.3 uses the following parame- 
ters: A = 540.0, B = 35.45, ki = -0.32354, a = 2.493, d = 1.420. 

are shown in Fig.|4] Spectral densities in the range 0.05 < k < 
0.3 h Mpc"' which is relevant for unbiased cosmological tests 
are plotted (dots) together with the theoretical wiggle function 
(Eqs. I1I7^ multiplied by an appropriate amplitude VA (solid 
line). We see that even for very large galaxy samples filling 
a box-Uke survey volume with a length of L = 500 /i"' Mpc 



8 



Ralf S. Koehler et al.: Probing Dark Energy with Baryonic Acoustic Oscillations 



error bars are mainly determined by sample variance and are of spectrum with a spectral index n — I (solid line), an index of 



the same order as the amplitudes of the BAOs. Therefore, only 
surveys with significantly larger survey volumes can show the 
BAOs with higher significance. 







Fig. 4. Extracted baryonic oscillations of dark matter particles (dots) 
with error bars including shot noise and sample variance at different 
redshifts ranging from z = 1.8 (upper left panel) to z = 3.8 (lower 
right panel). Superimposed is the theoretical wiggle function (Eq.0 
with an appropriate amplitude. 



We further see that the amplitudes of the wiggle functions 
in Fig.|4]decrease with increasing redshift, due to growth sup- 
pression (we verified the amplitude change with z according to 
standard linear structure growth). This is a direct consequence 
of the fact, that distortions of the power spectrum are directly 
translated into the boosted transfer function and thus into dis- 
tortions of the amplitude of the BAOs. To see how this works 
and how deviations from the assumed primordial spectrum with 
n - I affect the wiggle function, we rewrite the primordial 
spectrum in the form A:'^^'^ with |e| small compared to unity. 
For the assumed power spectrum, P(k, z, b) - A{k, z, b)k"T^{k), 
this choice gives 



J^^^-F(k,b,z)^ yfMk:Zb)k^^T^{k), 

V k Qni 



(8) 



where the phenomenological function (Eq.|6j effectively fits the 
amplitude and shape-dependencies (see Eqs.|3l3 



F(k,z,b)^ ^A{k,b,z)k' 



Lie iih ~ 

-^T,{k) + -^n(k) 



(9) 



It is seen that the phases of the wiggle functions remain un- 
changed by the FITEX procedure. This result is shown in Fig. 
|5] where the extracted oscillations from power spectra with dif- 
ferent spectral index models in the primordial power spectra are 
plotted. The oscillations were retrieved by calculating a power 



n - 0.9 (long dashed line) and a simple running spectral index 
model Hn*^ (short dashed line) and applying the FITDEX 
method. The parameters are chosen to conservatively cover the 
range of primordial power spectra not excluded by WMAP 3rd 
year data (see Bridges et al. 2006). 




0.1 0.2 

k [h-i Mpc]-i 



0.3 



Fig. 5. Wiggle extraction from power spectra with different primor- 
dial k dependencies. The resulting oscillations using FITDEX on pri- 
mordial models with a spectral index n = I (solid line) is compared to 
a model with different spectral index n = 0.9 (long dashed line) and a 
running spectral index with a = 0.1 (short dashed line). 



Only the amplitudes of the wiggle functions are modu- 
lated by possibly scale-dependent factors. Although a lot of 
interesting information about galaxy formation and cosmology 
is contained in the amplitude factors of the wiggle function, 
we skip this detail here and perform an amplitude rescaling, 
y/A(k,z,b)k'Q.^/a^ ^IA{k,b,z), with A(k,z,b) being a ba- 
sically arbitrary function after rescaling, summarizing distor- 
tions like structure growth, galaxy biasing, redshift space ef- 
fects, deviations from the Harrison-Zel'dovich case etc. Thus, 
with FITEX we are not deriving a pure undistorted wiggle func- 
tion, but a spectrum which is subject to all the distortions that 
work on the power spectrum. 

While the oscillations extracted from the Hubble Volume 
Simulations are consistent with the theoretical template, the er- 
ror bars prevent a precise evaluation of the accuracy of FITEX 
(all 4 cubes will be combined later, see Fig. ll5t . To get a bet- 
ter estimate of non-linearities that might wash out wiggles at 
higher k values, we use a non-linear analytic form of the power 
spectrum. This was computed by Jeong & Komatsu (2006) us- 
ing 3rd-order perturbation theory, which is able to model non- 
linear N-body matter power spectra to better than 1 % at z > 2 
mdk< 0.3/iMpc"'. 
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Fig. 6. BAO extracted from an analytic power spectrum of Jeong & 
Komatsu (2006) that includes non-linear growth at z = 2.3 (dots). The 
continuous line is the theoretical wiggle function derived in Sect.|2| 
multiplied by an appropriate amplitude. See the main text for a more 
detailed discussion of the treatment of the amplitudes of the wiggle 
function and the analytic power spectrum. 

Figure |6] shows the results of the FITEX extraction ap- 
phed to the analytic non-linear power spectrum using Hubble 
Volume parameters. As expected, the amplitude deviates from 
the theoretical prediction at small A:-values, where the shape 
contribution of the baryons dominates (see. Fig. but the 
phase of the oscillations in the analytical non-linear function is 
recovered to better than 1 % accuracy. This clearly shows, that, 
at least in theory, FITEX is able to model non-linear effects and 
recover the oscillations with high accuracy. 

As an additional test, the power spectrum of the Hubble 
Volume Simulation (OCTANT data), a wedge covering on 
eighth of the sky from z - 0-1.3 including redshift-space 
effects, was computed with about 60 million particles in a fast 
Fourier box with a legnth of 1600/? ' Mpc and the FITDEX 
algorithm used to extract the oscillations (Fig.^}. The oscilla- 
tions can be seen, although non-linear growth is heavily wash- 
ing out the wiggles at ^ > 0.22 h Mpc"' . 

Furthermore, we see a very good correspondence between 
a FITEX-extracted wiggle function obtained from N-body sim- 
ulations as well as the analytic non-linear power spectrum 
and the theoretical wiggle function of FITEX (Eq.^ in the 
HETDEX redshift range. Both linear and non-linear structure 
growth change the power spectrum and thus the BAOs, but 
FITEX still gives quite stable and easy to analyze results for the 
2nd to 5th BAO. The presence of the BAOs is consistent with 
the results of recent simulations of e.g. Springel et al. (2005) 
and Seo & Eisenstein (2005). However, no specific method to 
extract the BAOs was applied and their simulations did not in- 
clude lightcone effects. In the following we concentrate on the 



Fig. 7. Extracted BAOs using the Hubble Volume Simulation cover- 
ing one eighth of the sky from redshift to 1.3. Error bars include 
sample variance and shotnoise but do not cover washing out effects 
introduced by the correlation of modes through non-linear growth and 
distant observer approximation. 

effects of redshift space and scale-dependent galaxy biasing on 
the results obtained with FITEX. 

5. Extracting BAOs in Redshift Space 

The main problem in relating redshift space coordinates 
with comoving coordinates is the differentiation between red- 
shifting due to the expansion of the Universe and red-shifting 
due to peculiar velocities of the measured particles. Two effects 
might be important. 

(1) In the linear regime and in combination with the dis- 
tant observer approximation, peculiar velocities introduce a 
boosting factor p along the LOS. Only the A:-modes parallel 
to the LOS, k\\, are boosted by the square of the factor 1 -H 
I3{z,b)kl I {kl + k]) With I3(z,b) = -d\nD{z)l[d\n{\ +z)b(k,z)], 
D{z) the linear structure growth, and b(k, z) the biasing param- 
eter (Kaiser 1987). 

(2) In the non-linear regime, velocity dispersion of galaxies 
in (partially) virialized structures cause the well-known 'fin- 
gers of God' effect in redshift space. We assume that the ef- 
fect is mainly restricted to small scales and is describable by a 
simple exponential damping of the power spectrum by the fac- 
tor exp(-Q'A:^). Here, the redshift-dependent parameter a gives 
a measure of the average pairwise velocity dispersion in the 
given redshift shell. 

The combination of (1) and (2) constitute one version of the 
standard dispersion model of peculiar velocities which relates 
the one-dimensional redshift and real space power spectra by 

r' 2 

P,(k,z,b) ^ P(k,z,b) dfi[l + P(z,b)^^] e-"^'"''>\ (10) 
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Fig. 8. Redshift space spectra divided by their corresponding real 
space spectra at different redshifts (filled areas). Overplotted lines rep- 
resent the best fit of the kernel of Eq. <10t . 



where yu is the cosine between the direction of the actual 
wavevector and the LOS, i.e., jik - k\\, and - k^^ + k\ (for 
more detailed models see Scoccimarro 2004). On large scales, 
the real space power and thus the wiggle function is boosted 
in a ^-independent manner (Kaiser factor). On small scales, the 
power spectrum is dampened in a A;-dependent manner by the 
factor 0.5 yfn trfiak) I ak, with erf the standard error function 
of statistics. In Schuecker, Ott & Seitter (1996) it is shown that 
Gaussian redshift errors can be described with the same model 
(for /? = 0). Eq. ( I10> suggests that only the amplitudes but not 
the phases of the BAOs can be changed by this type of redshift 
space distortion. 

Figure|8]illustrates the combined effect in more detail. Two 
redshift space power spectra, as obtained from the Hubble 
Volume Simulation, were divided by their corresponding real 
space spectra at different redshifts and plotted as filled areas. 
The kernel of Eq. ( I10> was fitted to this ratio and the best fit 
results superposed (black solid hnes). We found that the model 
is accurate on the 5% -level over the complete redshift range. 
For z > 2, the parameter a is virtually zero because no virial- 
ization had time to take place on the scales shown. As redshift 
decreases, galaxy clusters start to form and virialize and the a 
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Fig. 9. Ratios of boosted redshift and real-space transfer functions 
for different redshifts. As expected, the redshift space effects increase 
with decreasing redshift in a smooth manner. Increasing ratios at high 
redshifts are an artifact caused by the glass-like initial load of the N- 
body simulations. To illustrate the relation between redshift space ef- 
fects and amplitudes of BAOs, the ratio between BAOs and the trans- 
fer function is plotted in the upper panel. 

parameter starts to distort not the amplitude in the two dimen- 
sional power spectrum depending on the angular component, 
like p does, but changes the amplitude depending on the mod- 
ulus of the k vector This means that the shape of the wiggles 
could be influenced by this effect as the amplitude of the power 
spectrum and thus the amplitude of the wiggles is dampened at 
high A:- values. 

Fortunately, this anisotropy only gradually changes the am- 
plitude of each consecutive mode over which is averaged in 
a ^-shell from perpendicular to parallel modes. This is illus- 
trated in Fig. |5] where the ratios of boosted transfer functions 
are plotted. A rough idea about the relation between redshift 
space effects and the amplitudes of BAOs is provided by the 
ratio between BAO and the transfer function plotted in the up- 
per part of Fig.|9] 

Two examples of wiggle functions boosted by redshift- 
space effects and extracted with FITEX are plotted in Fig.llOl 
The result of the anisotropic redshift space is that modes 
with a more parallel component get more weight than those 
with a more perpendicular component, which are not boosted. 
Whereas the overall shape of the wiggle function appears to be 
slightly changed, the phases of the wiggles are not distorted at 
all. Instead, the whole wiggle function receives a boost of the 
monopole component of the Kaiser factor. 

The fitting function is obviously able to model the redshift 
space distortions on the sub-percent level and oscillations can 
be derived. If the wiggles were not extracted with FITEX-like 
methods but by dividing the observed power spectrum by a 



Ralf S. Koehler et al.: Probing Dark Energy with Baryonic Acoustic Oscillations 



11 



i 




0.2 

k [h-i Mpc]-' 



0.3 



i 




0.2 

k [h-i Mpc]-i 



0.3 



Fig. 10. Comparison of wiggle functions at z = 1.2 (upper panel) 
and z = 2.3 (lower panel) obtained from real space (filled circles and 
continuous lines) and redshift space (open circles and dashed lines). 
Note the overall boost of the wiggle functions by redshift space effects 
which, however, leaves the phases of the BAOs invariant. The non- 
linear growth effects, start to wash out the wiggles at ^: > 0.2 h Mpc"' 
in the upper panel. 



standard non-oscillatory linear theory power spectrum, as was 
done by in Springel et al. (2005, Fig. 10) and Cole et al. (2005, 
Fig. 13), the baseline of the oscillations would start to depart 
from zero at smaller scales. Such an oscillation spectrum could 
not easily be used in a cosmological test as it fit only nontrivial 
theoretical model spectra. 

6. Extracting BAOs from Biased Samples 

In the Hubble Volume, CDM particles are simulated. What 
is actually observed are galaxies made out of dark matter 



and baryonic matter, especially in the case of the proposed 
HETDEX survey, where Lyman-a emitting galaxies are used 
as tracers of the underlying matter distribution. Because of the 
diff'erent properties of dark matter and baryonic matter (Benson 
et al. 2000, Le Delliou et al. 2006), both need not necessarily 
follow exactly the same density distribution. 

To test the stability of FITEX under such conditions, we 
introduce some sort of biasing into the N-body simulations. 
Formally, biasing is described by the mean biasing function 
(Dekel & Lahav 1999), 



b{6d)6d = {6g\6d) 



d6gP{6g\6d)Sg 



(11) 



that relates the dark matter distribution 6d to the galaxy dis- 
tribution Sg (Somerville et al. 2001). Here, P{Sg\ 6d) is the lo- 
cal conditional biasing distribution, i.e. the probability that a 
certain matter density contrast 6d with a variance cr^ is con- 
verted into a galaxy contrast 6g with a variance cr^. The bi- 
asing function can also be expanded into the moments, b - 
{b(6d)6l}/cr^, and P - {b^(6d)Sl}/(rl, with the biasing scatter 
o'l = - {6g\6d)f')lo-],- local, non-linear and stochastic 
biasing relation can be described by these moments to second 
order While b/b is a measure for the non-linearity and o-i,/b a 
measure of the scatter, the parameter b characterizes the linear 
biasing such, that if b - b, then Sg(x) = bSdix) (Wild et al. 
2005). 
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Fig. 11. Ratios of biased and unbiased boosted transfer functions for 
two samples with different scale-dependent biasing at z = 2.3. A factor 
of 2 was subtracted from the strong bias (dashed line) for enhanced 
visibility, while the weak biasing (solid line) remains unchanged. The 
theoretical wiggle function is plotted on top of the figure to illustrate 
the scale of the biasing effects. 



We use the simple stochastic and slightly non-linear 
Eulerian biasing scheme (e.g. Cole et al. 1998, Yoshida et al. 
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Fig. 12. Boosted wiggle functions extracted with FITEX and theo- 
retical wiggle functions for the unbiased sample (filled circles and 
continuous line), for the strongly biased sample (triangles and shot- 
dashed line) and for the less biased sample (open circles and long- 
dashed line). The theoretical wiggle functions are multiplied by the 
fo-moment of the biasing function. 

2001). First, local overdensities for each dark matter particle 
are computed by smoothing the distribution with a Gaussian, 
expi-r^/lr^) and ro = 3 h^^ Mpc. The overdensities 6d are 
then transformed into normalized overdensities, v - ddlcrd, 
with (Td being the root mean square of the ^^/-values. Finally, 
a function is introduced that describes the probability of a dark 
matter particle with overdensity 5d to be tagged as a 'galaxy', 
p{y) ~ (y — VcY, with Vc being a threshold at which the prob- 
ability is set to zero. This simulates the lower probability of 
galaxies forming in low density regions like voids. The param- 
eter a controls the stochastic spread of the biasing function as 
well as its non-linearity. A value of ff = 0.2 is a good com- 
promise between reasonable non-linearity and a comparably 
tight spread. Two samples with Vc - -0.2 and 0.6 were gen- 
erated. Both parameters, a and Vc, were adjusted to model the 
biasing parameters and number densities found by observation. 
SUBARU data (Hamana et al. 2004) suggests a comoving LAE 
density of « = 0.00301 ± 0.00194 Mpc"'' and a biasing pa- 
rameter ofb- 4-5 at a redshift z - 5, which should decrease 
slightly with redshift. As discussed in Sect.l^the biasing func- 
tion b{5d) should translate into a boost of the power spectrum 
b^(k,z)- Because the average slope of the biasing function is b, 
to first order we expect Pbias(fc) = P{k)b^{k,z) ~ P{k)P(z). 

The two prescriptions lead to a biasing that either decreases 
or increases with k (see Fig.ll 1>. Both values are consistent 
with the calculated b parameter. The comparatively strong scale 
dependence of the biasing becomes evident when the ratios of 
the boosted transfer functions are compared with the ratio of 
the wiggle amplitude to the transfer function plotted on top of 
Fig.[II] 



Figure^]illustrates, that FITEX extracts BAOs even from 
biased galaxy distributions. Plotted are the wiggle functions 
of the less strongly biased sample (open circles), the more 
strongly biased sample (open triangles) and the unbiased sam- 
ple (filled circles) at z = 2.3. The unbiased and biased oscilla- 
tions vary because of shot noise, as the biased samples have a 
2 and 10 times lower point density, respectively. Superimposed 
are the theoretical wiggle functions (Eq.Q, normalized to the 
unbiased wiggle function (continuous line) and multiplied by 
the first moment of the biasing, b (dashed lines). The boost 
of the power spectrum obviously translates into an amplitude 
boost of the BAOs. As in redshift space, only the phenomeno- 
logical fitting function as used in FITEX extracts the BAOs 
with a zero baseline, unlike when using a linear theory power 
spectrum. 

However, not only the overall boost of the power spec- 
trum, but also the ^-dependency of the biasing function should 
directly affect the amplitude of the BAOs. We tried to ver- 
ify this by fitting the Tu,{k) function supplemented by a 'lin- 
ear biasing model' to describe the extracted wiggle function 
Tsx{k) = (fl() + fli k)T„{k). We notice that both parameters aq 
and fli are highly degenerate, mainly due to the fact that only 
data points at the extrema of the wiggle function significantly 
contribute to the fit. Data points at or near the nodes of the wig- 
gle function are very insensitive to changes of the amplitude. 
Due to the degeneracy, any amplitude aq can easily be com- 
pensated by an appropriate value of the parameter oi. Though 
the full fits had the expected slopes, the effects were very small 
and of no significance for the cosmological test. We may thus 
conclude that in a cosmological test, which is mostly about de- 
tecting phase-shifts, it seems to be enough to fit the theoretical 
wiggle function with only an amplitude parameter, when the k- 
dependent biasing is as weak as simulated here. However, more 
realistic simulations are needed to analyze higher order effects. 

7. Cosmological Tests with BAOs 

The cosmological test compares the phases of the extracted 
BAOs with the theoretical wiggle function, projected onto the 
hypersphere of the observer As the theoretical template does 
not model the amplitude of the oscillations with high accuracy, 
we are not able to include amplitude effects in the cosmologi- 
cal test and marginalize over the wiggle amplitude. Future tests 
should include this parameter, which is sensitive to w by e.g. 
structure growth. 

A general feature of the test is that one first has to assume a 
certain reference cosmology to get the metric scale of the wig- 
gle function from the data, and to perform in a second step the 
comparison of the metric scale of the theoretical wiggle func- 
tion relative to the chosen reference cosmology (Glazebrook 
& Blake 2005). In the present paper, we are interested in the 
performance of the BAOs extracted with FITEX in constrain- 
ing the redshift-independent part of the w parameter of the dark 
energy. The amplitude of the theoretical function is also fitted 
but later marginalized. The remaining cosmological parameters 
are fixed to their values given at the end of Sect.^ 

As metric scales, comoving distances of the BAOs par- 
allel and perpendicular to the LOS are used. Recall that in 
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the plane-parallel or distant observer approximation, comov- 
ing distances between two points with the same angular posi- 
tions but different redshifts are jcy = J]~" ^dz - -gpy, 
where H{z) = HqE(z) is the Hubble parameter at redshift z and 
E{z) the w-dependent transformation of the Hubble constant, 
Hq, from redshift zero to z (e.g., Peebles 1993). To compute not 
distances between two points in space with different redshifts, 
but parallel scaling factors between the two sets of cosmologies 
denoted by primed and unprimed symbols, one has to take the 
limit Az-Z2-Z\ — » 0, which yields 

dx! 

x[, =X|| — . (12) 

Comoving perpendicular separations represent distances two 
points with the same redshift but different angular positions are 
separated, x_i_ - j; dz, where 6 is the angle between the 
two points. The scaling relation for the perpendicular compo- 
nent between two cosmologies thus becomes 

, x' 
x^ ^ xj_ — . (13) 

X 

In Fig. [0]the transformation ratio of the wiggle function is 
plotted at different redshifts for two cosmologies that differ 
from the reference model only in vvq (the redshift-independent 
part of w). The reference model is the concordance cosmol- 
ogy. Fig.[Ol shows that the largest deviations from the refer- 
ence model appear at redshifts in the range 0.5 < z < 1.5. 
Differences in parallel and perpendicular scaling introduce 
anisotropics that could be measured, too (see below). 
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Fig. 13. Transformation of perpendicular (solid) and parallel (dashed) 
distances in comparison to the concordance model with wq = -0.7 
(upper part) and Wo = -1.3 (lower part). 

FigurelT^ shows the scaled wiggle function at redshift z - 
2.5 for Wo = -0.7. The reference and the distorted wiggle func- 
tion are plotted. Note that an elongation of the coordinates in 
real space leads to a compression of the coordinates in ^-space 
and vice versa. In the upper panel, it can be seen that the wig- 
gle function gets more and more out of phase as k increases. 
Thus the sensitivity of a cosmological test should increase at 
high values of k. However, the amplitudes of the oscillations 
decrease exponentially with k (e.g.. Silk damping). This can 
be seen in the lower panel of Fig. where the scaled wig- 
gle function is subtracted from the reference (wq = -1) wig- 
gle function: Though the difference in phase is increasing, the 
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Fig. 14. Scaled wiggle functions at z = 2.5 for vvo = -0.7 (dashed 
line) compared to the wo = -1 cosmology (solid line, top panel). 
The normahzed modulus of the difference between the reference cos- 
mology and the scaled cosmology is shown in the bottom panel. The 
dashed vertical hues mark the region relevant for the cosmological 
test. 



amplitude decreases so fast that the difference between the pro- 
jected oscillations increases at small values of k and decreases 
at high values of k. 

The present version of the cosmological test comprises the 
following steps: (1) Angular and redshift coordinates of the 
'observed' (simulated) data are transformed into comoving co- 
ordinates using a concordance model as the reference cosmol- 
ogy. (2) The comoving data sample is Fourier-transformed to 
estimate the boosted transfer functions. (3) BAOs are extracted 
with FITEX. (4) The theoretical wiggle function is computed 
for the corresponding cosmological model. (5) The parallel and 
perpendicular scaling factors between the test cosmology and 
the reference cosmology are computed at the redshift of the 
data sample. (6) The theoretical wiggle function is scaled by 
constructing a three-dimensional wiggle function, scaling each 
mode and collapsing it to one dimension. (7) The observed 
wiggle function is compared to the re-scaled theoretical wig- 
gle function assuming Gaussian statistical errors {x^ test). Only 
the steps (4-7) have to be repeated to test a different cosmology, 
allowing a fast test of various parameters. 

Note that the three-dimensional wiggle function is only 
constructed to simulate the effects of different parallel and per- 
pendicular scaling factors. It does not address the anisotropy 
introduced by redshift space effects, because we found that this 
only boosts the various modes but not change the phases of the 
theoretical wiggle function (see Sect|5}. The boost effect of the 
redshift space distortions, along with other boosting effects like 
biasing and growth suppression, will force us to add an addi- 
tional parameter to the cosmological test. 
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Figure[2]shows in the upper left panel the combined wig- 
gle function of the four XW wedge cubes of the Hubble Volume 
Simulation. The combination is performed by averaging over 
^-bins of discrete size. The errors represent the standard de- 
viation of all wiggle points in the corresponding ^-bins. After 
averaging, the 2nd to 5th BAOs are now clearly visible. We thus 
conclude that a volume at least 4 times larger than SDSS has 
to be sampled in order to see the BAOs with FITEX with high 
significance. One has to keep in mind, though, that our cubes 
are located at much higher redshifts, where the smoothing ef- 
fect of non-linear structure growth is obviously much smaller. 
We also neglected the fact, that larger sky coverages reduce the 
sizes of the fundamental modes, leading to higher (better) sam- 
phng rates of BAOs. 
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Fig. 15. Combined wiggle functions from the four Hubble Volume 
cubes. The upper left panel shows the combined data point (dots) with 
self-consistent error bars for the 100 square degree field. The upper 
right panel shows expected error bars for 200 square degrees and the 
lower row with error bars for 300 and 600 square degrees, respec- 
tively. The continuous lines are the theoretical wiggle functions for 
the reference (concordance) cosmology. 



With the 100 square degree sample, the possibility of data 
points appearing at positions which severely inhibit the accu- 
racy of the cosmological test is relatively high, due to the high 
sample variance. Especially as only a few data points, which 
are far away from positions where theoretical wiggle functions 
of different cosmological models intersect, have great statistical 
weight in the cosmological test. This fact is shown in Fig. 1141 
where the modulus of the difference between the theoretical 
transfer functions for two models depending on position in k- 
space is plotted. 

Figure ^] shows the results of the cosmological test ob- 
tained with BAOs extracted with FITEX under different sur- 
vey conditions. For surveys covering 100 square degrees we 
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Fig. 16. distributions for the wq parameter obtained after marginal- 
ization over the amplitude parameter and for different conditions of 
BAO extraction. Short-dashed line: dark matter distribution in config- 
uration space over 100 sq. deg obtained from the simulation of dark 
matter particles in configuration space. Continuous line: biased parti- 
cle distribution in redshift space over 100 sq. deg. Dotted line: biased 
particle distribution in redshift space over 300 sq. deg (extrapolation). 
Dashed-dotted line: biased particle distribution in redshift space over 
600 sq. deg (extrapolation). The horizontal dashed line marks the Icr 
error level. 



would expect eiTors of cr„,^ - 0.25. The asymmetry of the error 
distributions directly reflects the increased sensitivity of metric 
scales on w for wq > -1 compared to wo < -1. At these wq 
values, the theoretical wiggle function gets quite compressed, 
so that at large k values the amplitude decreases rapidly due to 
Silk damping. When wq — > -oo, the theoretical wiggle func- 
tion formally mimics the case of zero amplitude, so that low 
values of wq are excluded with about the same significance as 
the "no-wiggle" case of zero amplitude. The minima of all er- 
ror distributions plotted in Fig.^]are slightly off-set from the 
input W() = -1 value, but well within all formal Icr ranges. 

Note that the constraints given above are based solely on 
the phases of the oscillations in the one-dimensional power 
spectrum, assuming perfect knowledge of the sound horizon at 
drag epoch. The w-constraints are expected to further improve 
when the full two dimensional redshift space effects as well as 
the amplitudes of the BAOs are used. 

Perhaps the most important result of the present investi- 
gation is the comparatively small decrease of only 8% of cTh,^ 
when comparing the w test with unbiased matter distributions 
in configuration space, (dashed line in Fig.ll6> i.e., the simplest 
survey condition, with the test that includes all aforementioned 
effects (solid line). The robustness of the FITEX method 
to extract BAOs from boosted transfer functions with phe- 
nomenological fitting functions is mainly related to the fact 
that only the amplitudes of the wiggle functions are boosted 
by redshift-space and biasing effects, but the phases are not 
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affected. 

To give an idea of the performance of FITEX based on sur- 
veys of larger size, we extrapolated the results from the 100 
square degree data. The upper right panel of Fig.^jshows wig- 
gle functions extracted from the 100 square degree data but 
with error bars expected for 200 - 600 square degrees. The er- 
ror bars were extrapolated by scaling them according to the 
expected number of additional fc-modes. The different survey 
fields are thus assumed to be not adjacent. 

Figure [21 shows that increasing the survey volume to 200 
square degrees does not reduce the error bars significantly, still 
allowing for crucial data points at peaks of the oscillating wig- 
gle function to be reduced to zero amplitude. Error bars are still 
of the same order of magnitude as the oscillations. Only at 300 
square degrees (lower left panel), the probability for smoothing 
out the oscillations due to statistical errors gets considerably 
low. Following our extrapolations, for a survey of 600 square 
degrees, the errors are expected to get much smaller than the 
amplitude of the oscillations and the chance of statistical devi- 
ations from the predicted function is of an order of magnitude 
that allows for very accurate cosmological probing. 

The accuracies of the wq parameter (Fig. ll6t range from 
(T„„ = 0.25 for 100 sq. deg which corresponds to 3 SDSS 
volumes, to cr^^ - 0.09 for 600 sq. deg, corresponding to 21 
SDSS volumes. These results are obtained for biased samples 
in redshift space, i.e., the most complex survey condition. The 
improvement is directly related to the larger number of inde- 
pendent ^-modes in larger survey volumes and is consistent 
with results obtained with Monte-Carlo simulations (Blake & 
Glazebrook 2003, their Fig. 8). 

Future studies should also investigate the various parame- 
ters of the cosmological test to get constraints under less strong 
priors. 

8. Discussion and Conclusions 

We have seen that it is critical for the extraction of BAOs to 
either fit a phenomenological function to the continuum of the 
boosted transfer function (FITEX), or to calculate the shape 
of the boosted transfer function, including all distortions down 
to the sub-percent level. The latter approach has been carried 
out recently. Glazebrook & Blake (2005), Hutsi (2005), and 
Springel et al. (2005) used a linear theory power spectrum 
and divided the 'observed' data by this reference. Almost the 
same approach, but with the addition of linearizing the observa- 
tional data before the division, to correct for non-linear struc- 
ture growth was used in Angulo et al. (2005). There are two 
drawbacks of this method: 

First, the ratio of a complex power spectrum and a refer- 
ence power spectrum does not resemble a wiggle function as 
defined in the text, but something more complicated (see Eq.|^. 
The simple form of the theoretical wiggle function (Eq.^ is 
strongly deformed by both the chosen transfer function of the 
reference spectrum and the distortions introduced by redshift 
space and galaxy biasing. These distortions make the compar- 
ison with theoretically expected model spectra unnecessarily 
complicated. FITEX transforms the observed power spectrum 



into a boosted transfer function where the subtraction of a non- 
oscillating phenomenological continuum function is directly 
related to r„,(fc), multiplied by an amplitude factor 

The second drawback is the linear theory power spectrum 
itself. By computing a correct shape of the power spectrum, 
one has to rely strongly on assumptions about the nature of the 
dark matter and various other physical processes. The method 
is no more the 'assumption free' approach that was intended. 
Furthermore, the calculation of the amplitude factor that in- 
cludes all the distortions ranging from linear structure growth, 
non-linear structure growth to redshift space distortions and bi- 
asing is a problem. This calculation has to be accurate to the 
sub-percent level, as the oscillations themselves make up only 
~ 2% of the transfer function. 

A more phenomenological approach was conducted by 
White (2005). He assumed that redshift space distortions as 
well as linear and non-linear structure growth could be calcu- 
lated with reasonable accuracy. Leaving only the ratio between 
dark matter power spectra and 'galaxy' power spectra, he found 
that one can fit this ratio with about 2% accuracy using a cubic 
polynomial. This adds further evidence that a phenomenolog- 
ical approach can model power spectrum distortions without 
assuming much about the physics of the distortions themselves. 

With the FITEX method, we have shown that one is able 
to extract the BAOs without any strong priors on the nature 
of dark matter or the shape of the power spectrum from the 
complex multi-component transfer function. All the aforemen- 
tioned effects are modeled by the continuum function and need 
not be known very accurately, as long as the phases of the os- 
cillations are not distorted or washed out. This extraction is ac- 
curate on the sub-percent level within fc-ranges relevant for the 
cosmological test. Furthermore, FITDEX is able to disentangle 
phase information from amplitude information which could in 
the future be used to test the geometry of the Universe on the 
one hand and the growth history on the other hand. This type 
of cosmological test could be used to discriminate Dark Energy 
from Modified Gravity theories. 
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